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(54) Improvements in or relating to signal processing 



(57) A method and apparatus for reducing noise 
correlation in a partial response channel through optimi- 
zation of a look-ahead maximum likelihood (ML) detec- 
tor. In the method of the present invention, the ML 
detector is optimized in light of the noise correlation 
generated by the partial response channel. The 
in^proved ML detector provides comparable perform- 
ance to. or better performance than, a Viterbi detector in 
the presence of colored noise. In the present invention, 
a set of finite impulse response (FIR) transversal filters 
are used as the ML estimator for the look-ahead detec- 
tor. The weighted sum outputs of the FIR fitters are com- 
pared to a set of thresholds based on previously 



detected data to make the decision for current detec- 
tion. The present invention improves the ML detector's 
performance and reduces its conrplexity by optimizing 
the coefficients of the FIR filters in the>presence of the 
correlated or colored noise. The SNR (signal-to-noise 
ratio) of each FIR filter is determined for a range of coef- 
ficients based on the noise autocorrelation of the chan- 
nel for a given user density, and the coefficients 
providing the highest SNR are selected for each deci- 
sion function. The result is a noise whitening ML detec- 
tor providing improved performance and lower 
complexity than prior art ML detectors. 
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Description r.r::.-. . 

FIELD OF THE INVENTIGM : = . .c-.\ — :y: It rv . - : ' ; ; : . - 

5 The present invention relates to the field of signal processing in a channel, and, in particular, to signal processing 

in a partial response channel; * ■ ■ ■ > . . r . • :. . . . 

BACKGROUND ' ^ / : : : : v ■ : ^ ^ 

10 In magnetic recording devices, such'as magnetic disks and tapes, a recording head is used to read and write infor- 
mation to arid from a magnetic surface. In a typicar rotating medium-based storage system,.data is stored on magnetic 
disks in a series of concentric Vacks". These tracks are accessed by a read/write head that detects variations In the 
magnetic orientation of the disk surface. The read/write head moves back and forth radially on the disk under control of; 
a head-positioning servo mechanism so that it can be selectively positioned over a selected one of the tracks. Once in 

75 position over a track, the servo mechanism causes the head to trace a path that follows the center line of the selected 
track. . ' : ■ . ^ • : ■ 

Generally, the inductive' recording head consists of a slit toroid made up of high permeability magnetic material and 
wound by several conductor turns. The toroid contains a gap which is positoned over the data tracks on the magnetic 
recording surface. To record, a current is generated through the conductor windings, altering the magnetic field in the 

20 toroid. At the Ideation of the gap, the amplitude of the magnetic field is targe enough to record on the magnetic material 
of the storage device to a sufficient depth. The amplitude of the magnetic field falls off sharply away from the gap. 

By manipulating the current through the conductor windings, the magnitude and direction of the magnetic flux at 
the location of the gap can be modulated in such a fashion as to encode Information into the magnetic surface of the 
storage device. A pattern of external and internal fields are created as the head and recording surface are moved rela- 
ys five to each other. These patterns are similar to a series of bar magnets of changing polarities. The polarity transitipns 
are then readable as transitions in the magnetic flux at the recording surface.* In read mode, as the magnetic storage 
surface moves across the gap in the head, the magnetic field of the storage surface is detected at the gap, and a voltage 
is induced in the coil proportional to the rate of change of the flux. The read channel then processes this analog voltage 
signal, typically converting the flux transitions into a series of pulses, and converts the analog signal into digital data. 

30 Magnetic storage devices sometimes use analog peak detection to process incoming read signals. However, as 
recording density, or "user density" (Du), increases the analog peak detection scheme becomes unreliable because of 
the large amount of inter-symbol interference (ISI) between adjacent pulses. Alternatively, a partial response maximum 
likelihood (PRML) channel can be used to increase the user density. However, this method requires very good equali- 
zation of the read signal to shape the read signal into the desired partial response (P R) target waveform. 

35 User density (Du) is a measure of how the pulse width relates to the data rate in the form of a ratio of the analog 
signal pulse-width at 50% magnitude over the length of the data period. A higher user density means that signal pulses 
resulting from read transitions are spread over a larger number of dock cycles, increasing the interaction or overlap of 
adjacent pulses, and resulting in greater ISI;^ ^ 

A partial response chainnel does not remove ISI from the channel signal. Instead, the signal is equalized to provide 

40 a target pulse waveform having predefined values at specific sampling instants. The effects of ISI on the target pulse 
waveform are predicted and incorporated into a pulse detection scheme which determines the occurrence of a signal 
pulse based on the sample history of the waveform and model waveform prediction. The history of the sample wave- 
form is matched to a set of model waveforms, and the most likely model waveform is used to decode the signal into 
binary data. Hence, it is often referred to as "maximum-likelihood" detection. 

45 Partial response channels are therefore more robust than analog peak detectors in the presence of ISI. However, 
as described below, a problem exists in PR channels due to the equalization required to generate the target pulse wave- 
form, as well as the process used to generate a signal metric for determining the most \\ke\y waveform. The PR equal- 
izer and ML estimator used to generate the signal metric each serve to correlate noise in the channel. As a result ML 
detector performance is impaired by a reduction in signal-to-noise gain. 

50 The most commonly used maximum likelihood (ML) detector is the Viterbi detector, which is optimized for detection 
of signals in the presence of white noise. It is, however, somewhat complicated to implement especially for higher 
orders of PR channels (e.g. EPR4, EEPR4) and very inefficient for parallel/pipeline irrplementation architectures. Fur- 
thermore, its performance degrades as the noise become correlated. Look-ahead ML detector hardware requirements 
are usually less than that for a Viterbi detector, and look-ahead detectors are more suitable for parallel/pipeline imple- 

55 mentation. However, the performance of look-ahead ML detectors is usually inferior to a Viterbi detector. 

Look-ahead detectors are described in A. M. Paters article "A New Digital Signal Processing Channel of Data Stor- 
age Products, "Digest of the Magnetic Recording Conference, June 1991 , pp.E6-E7: in A.M. Patel, et al., "Performance 
Data for a Six-Sample Look-Ahead (1.7) ML Detection Channel." IEEE Transactions on Magnetics vol.29, no. 6. 
pp.401 2-4014, November 1993; and in Patel's U.S. Pat: No. 4.945,538. A look^ahead detector is also described in U.S. 
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Pat. No. 5.31 1,178 to T.Pan and R.G.Yamasaki. / ■.■ : , 

The Patel and PanA'amasaW channels use partial response signaling to limit the number of isolated pulse non-zero 
samples. The Patel channel Is equalized to EPR4 with three non-zero sanples, and the Pan/Yamasaki channel is equal- 
ized to EEPR4 (or E^PM) where the four non-zero samples, a, p. p and a, are 1.3,3 and 1 . respectively. Both channels 
5 use the (1.7) RLL code to simplify signal detection. . ci^-c : rlv fo i:.:;t :^ : : ;f:;':':.f . 

Run-length limited (RLL) codes are useful because they place upper and lower bound^ .pn the nurnber of dock 
cycles occurring between transitions. The upper bound is very important because dock recovery is based on the occur- 
rence of these transitions. For example, a long train of zeros in a data sequence produces no transitions: and the clock 
recovery circuit has no input pulse with which to synchronize its tracking. In this situation, the data recovery timing might 
10 drift out of phase. RLL codes ensure that sufficient transitions:occur for the clock recpyery circuit to maintain the correct 
timing phase and frequency Also. by:maintaining the lower t)Ound on:the number of zeros/between consecutive ones, 
signal pulses in the read channel are separated; to reduce the-Jntersymbor interference: (isi); c interaction 
between adjacent signal pulses:. ' h " ■ • .c rr r . j;: }■ : .a::?:;.: . e v * a . * r- ':\ -v ; ^ : . . ; • . 

The (1 i7) RLL code is characterized by a minimum of one. "Q" and a .maximum of seven "O's** between consecutive 
75 "1's". In the modified non-return-to-zero (NRZl):format,'wheraeach "rJ"/is=r5presented by a trjansition. and each "O" is 
represented by the lack of a transition, the (1,7) RLL code is sufficient for clock recovery purposes. Further, by main- 
taining the minimum of one **0" between consecutive *'T's" . transitfons are separated so a&:^to^beJCiifferentiable from one 
another. ' - ■ ' c-" ■ y. Kr-r.::^vr: -;c-:o: rr'"" ■ . 

Figure 1 illustrates, with solid curve 1G0.'an isolated Lorentzlan-pulse equalized and sampled! for EE PR4. and a 
20 minimum distance error event shown as a dashed curve 1 01. rljie pulse is the result of applying a.positiye unit step tran- 
sition, e.g. a transition from a magnetic read head, to a {1-D)(1+D)^::chanriel (EEPB4);.:The circles on. solid waveform 
100' show the isolated pulse samples for EEPB4: Samples of minimum distance^ error event 101 are indicated by tdia- 
monds. Minimum distance error event 1 01 is equivalent to the isolated pulse (tOO) shifted by one cpde period in^ime, 

In most ML detectors.' including the Viterbrdetector, the square of the Euclidean distance between two sequences 
25 is often used as the signal metric. The Euclidean distance is defined as.the squareToot of the surn of the squared^error 
terms between two waveforms at each sampling instant. The squared Euclidean.distance:er:ror of a sample combiriation 
is'defined as: " ,v.. ■ r ■ ; ■ • - - 

30 E^=(yo^ya)^ + (yi:-.yb):^.+ (y2> . ' [1] 

where ya. yb, yc, yd and ye are the expected valu^ for themodel. wavefbrni;-yo/ yi,:y2. 73 -and 74 are cprisecutive sam- 
35 pie values; and dj is the sample error at y|. i^. \. ; . : . / v e ^u . : : :r ^. : . ^ t* 

In Figure 1. the isolated pulse waveform TOO and the minrmuntdtstance error event 1.01 have sample, errors of val- 
ues 0,1,2,0,2 and 1 for samples yMi. yo. y^. y2. Vz and y4, respectively, where yi^^^ is referred. to as a look-back sample.. 
The squared Euclidean distance between sampled waveform 100 and the minimum distance error event 101 is labeled 
^M\i/' c*MiN^ equal to the sum of the squares.of these values. which adds up to a value of 10. 
40 In the ML detector of U.S. Patent No.5,3 1 1 . 1 78. thresholds for each decision function are determined by minimizing 
the error between an expected sanrrple model waveform and:the nearest valid sample model waveform. For this reason, 
the thresholds for each function are dependent on- the presentstate. as the next valid states are determined by the 
present state. < . ; ' : m . ^ ? • 

In the calculation of the Euclidean distance for use as a signal metric, :f or a f irst.model wavefbrrn to be selected over 
45 a second model waveform, r . - - 3 :. * . . . . ; , - 1 

' /■ VE2-^<>j7' ■ ^ [2] 

50 where E^^ is the squared Euclidean distance between thef inst model waveform and the sampled signal, and is the 
squared Euclidean distance between the second model waveform and the sampled signal. The squared terms of the 
sampled signal are canceled out to yield a linear decision function. This will generate the following threshold decision: 



55 



{ya2-yai)yo + (yb2-ybi)yi> (/Cg-yciOya + (yd2-ydi)y3 + (y^a-ye i)y4 

<'(ya2-ya i)(ya2> ya.,) + L. + (ye2-Yei)(ye2 + yej)' ' ' 



where yo* Ya » ya. ys are the actual sample values; yai , ybi . yei . ydi and yei are the expected values for tiie first 
model waveform; and ya2. yb2. yc2. ydg and ye2 are the expected values.forthe second model waveform. . . 
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Tor the sample points at which the model waveforms are equal, the products including' the factor (yk2 - yki). where 
k is the associated sample point, are canceled. In the negative phase, i.e. for negative transitions, the thresholds are 
the same, but the inequalities are r6versed.and the decisldn.functions are multiplied byM ^ : 

The peak check decision functions from Figure 5A of U.S. Pat. No. 5.311.178 are shown below for a=:1 and p=3. 
5 Fw (Fz in the patent) is normaliz«i with a sign change. ^' ^ •^• '^ :: ; :: : 

i ■ ' - ]a} Fx = Fy= d.5y + '-^^ - . " " 

^0 b. Fw = 0.5yo + yr-;0.5y3 + 0.5y4- [5] 

yO, yl , y2, yS and y4 are sequential code rate sanple values. The linear equations Fx; Fy and Fw have sample val- 
ues that are multiplied by coefficients which maximize the detection threshold between all possible signal waveform pat- 
terns and their minimum distance erro^waVefdrms, For exiample. each decision function (Fx, Fy, etc.) can be 
15 represented by the following equation: 

K 

• -^^^ F(n)=2c^^ - [6] 

where C| are the coefficients chosen to maxihnize detection, and k is the number of samples used in the ML detection 
scheme. . * 

From equation [6], it can be shown that each decision equation has the same form as an FIR filter Therefore, the 

25 detected noise and tiie detector hblse bandwidth are affected by the decision function. Further, in the prior art. decision 
functions are designed under tiie assumption tiiat the noise at the sampled input of the detector is additive white Gaus- 
sian noise (AWGN). However, partial response channels comprise filter and equalizer components for shaping the input 
waveform to the target PR pulse. The PR equalizer, in particular, typically supplies high frequency boost to provide 
pulse-slimming of the input pulse. . - 

30 The filter and equalizer components in the PR channel shape the noise of the system in the same manner as the 
input pulse, thereby coloring, or correlating, the input noise according to the filter and equalizer frequency response 
functions. The target puise signal at the input of the detector thus contains a noise component that is no longer white 
in nature, but has an effective noise bandwidth. For this reason, performance of the decision functions within the detec- 
tor IS reduced due to the correlation of noise. A partial response channel is desired that achieves optimum decision per- 

35 fbrmance in the presence of colored noise. 

SUMMARY OF THE PRESENT invention: 

The present invention provides a hnethod and apparatus for reducing noise correlation in a partial response channel 
40 through optimization of a look-ahiead majdmurn likelihood (ML) detector; In the method of the present invention, the ML 
detector is optimized in light of the noise correlation generated by the partial response channel. TTie improved ML 
detector provides comparable performance to; or better pigrformance than, a Viterbl detector in the presefnce of colored 
noise. : . : " 

In the present invention, a set of finite impulse response (FIR) transversal filters are used as the ML estimator for 
45 the look-ahead detector The weighted sum outputs of the FIR filters are compared to a set of thresholds based on pre- 
viously detected data to make the decision fdrcurrent detection. The present irrvention improves the ML detector's per- 
formance and reduces its complexity by optimizing the coefficients of the FIR filters in the presence of the correlated or 
colored noise. For a range of coefficient values, the SNR (signal-to-noise ratio) of each FIR filter is determined based 
on the noise autocon'elation of the channel for a given user density, and the coefficients providing the highest SNR are 
so selected for each decision function. The result is a noise whitening ML detector providing improved performance with 
lower complexity than prior art ML detectors. 

BRIEF DESCRIPTION OF THE DRAWINGS - - - := 

55 The ppresent invention will now be further descnbed. by way of example, witii reference to be accompanying draw- 
ings in which:- 

Rgure 1 Is a waveform diagram of an isolated EEPR4 pulse! with a minimum distance error event; 
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Figure 2 is a waveformidiagrarn of a. sampled worst-case EEPR4.vyayeform with a niiriirnum dl.starice error event; 
Figure 3 is a trellis state diagram.for a (1 J) RLL coded"EEPR4 loofeah^ 



Figure 4 is a simplified block diagram of a partial response (Pp) detection chanriel;,. 



Figure 5 is a decision table for positive and negative phase EEPR4 detection in one embodiment of a five sample 
look-ahead detector; 

10 Figure 6 is a state table for a ten state look-ahead ML detector; - - - - 

Figure 7 is a signal-to-noise, gain- table fpr yarious ML detectors: atyarying user densities;. - - 
Figure 8 Is a block diagram of a ML detector comprising FIR transversa! fi^^^ . ; . : , 

Figure 9 is a now diagram of a method for optimizing a estimator in the presence of colored noise; and 
Figure 10a frequency response diagram illustrating the correlating properties -of a partial response channel. 

20 DETAILED DESCRIPTION OF THE PRESENT INVENTION 

A method and apparatus for reducing noise correlation in a partial response channel is described. In the following . 
description, numerous specific details, such as EEPR4 channel detection and the variance of the noise power spec- 
trum, are described in order to provide a more thorough description of the inventioa It will be apparent; however. to,one 

25 skilled in the art, that the present invention may be practised without these specific details. In other instances,-wellr 
known features have not been described Jn detail in order not to unnecessarily obscure :the present invention. -.^ . ^ 

Prior art PR channels require complex look^ahead detectors to approach the signal-tO:noise gain performance of 
Vlterbi detectors under similar conditions.: Thus, in prior art circuits, look-ahead detectors with fewer look-ahead sam- 
ples are implemented with an associated reduction in performance. Furtiier, prior art ML detectors are prediced on the 

30 aSiSumption that noise at the input pf the ML detector is additive white Gaussian, noise. :r 

The method and apparatus of the present invention implement a lopk-ah^ad detector of fewer look-ahead samples - 
(e.g. four to six samples) in which the decision functions are optimized for the presence of polored noise. The expected 
correlation of the noise at the input of the ML detector is deterrnlned; and the; ML detector is designed to rnaximize the 
signal-to-noise gain of each ML decision function in light of the noise variance at the ML estimator output. Performance • 

35 approaching, and in some cases exceeding, that of a Viterbi detector of the prior art. is achieved. . ; 

Figure 2 shows a "worst case" signal waveform 200 with the associated minimum distance error event indicated by 
dashed curve 201. As shown in Figure 1, dMiN^ (minimum squared distance); for an isolated EEPR4 pulse is equal to 
ten. However, ten samples (yo-yg) are required to achieve a reading of dMiN^=10 for the worst case waveform of Figure 
2 due to the effects of ISI on non-isolated pulses. Thus, assuming AWGN (i.e. uncorrelated noise), a ten sample look- 

40 ahead detector is required.to achieve the maximum detector signal.-to-noise.(S/N) gain, .which is accomplished by max- 
imizing the minimum distance between waveforrns. . : 

A ten sample look-ahead detector is much more complex than ltie prior art four to six sanriple Ippkahead detectors. 
However for practical user densities (Du). the noise at the detector input Is not AWGN, but Is correlated (colored) by the 
frequency boost equalization required to achieve the desired partial response ^PR) target pulse signal. Correlation of 

45 noise in the channel acts to reduce the actual signalrto-noise gain of detectors operating under the .assumption of 
AWGN: Greater actual signal-to-noise gain is therefore possible with a less cornplexjook-ahead detector when the 
detector accounts for noise correlation. , ■: . ^^ 

The effects of the PR equalizer on the additive channel noise is best illustrated in the frequency domain where tiie 
effects on the noise bandwidth can be clearly seen.. Figure lOJs.a graph of the frequency responses for a seven-pole 

so equalizer function, a five-sample decision, function (FIR filter), and the combined, response for the series connection of 
the equalizer function and the decision function. TTie horizontal axis represents frequency in tens of megahertz. The . 
vertical axis represents gain. 

Curve 1000 is the frequency response for a seven-pole, two-zero filter representing the equalizer for a PR channel. 
The form of curve 1000 is that of a towpass filter with high frequency boost Noise in the PR channel is colored or cor- 

55 related by this frequency response so that tiie noise power spectrunri takes the shape of the equalizer response where 
the gain magnitude is squared. Curve 1000 is therefore representative of the noise spectrum at the output of the PR 
equalizer. 

Curve 1001 is tiie frequency response for a ML estimator decision function having coefficient values of [0.5.1 .0,-.1 
0.5]. Each decision function has an associated frequency response determined by its coefficients. The form of curve 
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1001 for this example is that of a bandpass filter dehtered near 7.5 MHz with a reduction in gain in the frequency region 
cbrresponding to the high frequency boost of cur/e 1000. Any signal processed by the^ML'esti decision function 
is shaped by this^ fi'equericy rissponse/ including- noise- lTie1o^ of curve 1001 cah theriefore be used as a whitening 
filter with respect to colored noise at the input of the ML e^^ 

£ Curve 1002 is the combined frequency response generated by the PR equalizer aind ML estimator in series. The 

shape of the curve Is deterniined by the product of the maghittides of curves 1000 and 1001 at each frequency." The 
noise power spectrum for- the P R channel including the PR equalizer and the ML estimator is therefore shaped by the 
square of the magnitude of curve 1002. For this example, the combined frequency response is a bandpass function with 
a center frequency near 10 MHz. By altering c6eff^eients of the ML estimator, the effective noise bandwidth associated 

10 with the combined frequency response can be i-edutied. Thereforei the present invention optimizes the coefficients of 
the ML estimator to rediice the effective noise bandw at the ML estimator output, thus whitening or decpn-elating the 
correlated noise component. As the rriinimtirh distance error is also dependent on the ML estimator coefficients, decor- 
relation of noise must be balanced with rfiaint^hahee of the minimum distance error to achieve the highest slgnal-to- 
nolse gain, 

75 In the following description, a (1 .7) RLL code EEPR4 channel ideally equalized for a series of Lorerrtzian pulses is 
used as an example. The invention is desaibkl iri relation to a (1;7) EEPR4 look-ahead detector embodiment. It will be 
apparent that the present invention can be' applied to other parlial resiDOhise channels as well. Further, a trade-off exists 
between detector performance and complexity. So as not to obscure the invention, a simpler five sample look-ahead 
detector will be described. 

20 Rgure 4 shows a block diagram of the PR channel: Automatic gain control (AGC) amplifier 401 sets the input signal 
amplitude so that input signal 400- is samplied at preset levels. AGC amplifier 401 provides the gain control signal 402 
to PR equalizer 403. Signal 402 comprises Lorentzian pulses generated from signal transitions, such as those in a mag- 
netic read channel. PR equalizer 403 equalizes the signal to the desired partial response waveform, performing pulse 
slimming on the input waveform by providing high frequency boost. Ideally, equalized signal 405 generated by PR 

25 equalizer 403 contains pulses whose sample! values conform to expected waveform values. Signal 405 is provided to 
maximum-likelihood (ML) estimator 407 within ML detector 406. - 

Within ML detector 406. ML estimator 407 processes the sample values through a set of decision functions to gen- 
erate a set of detection values, or metrics, labeled as signal 40a These detection values are provided to decoder 409 
for comparison with state dependent detection threshold values. Decoder '409 determines the detected data value and 

30 provides the' detected data as output signal 410. Concurrently, decoder 409 provides update signal 41 1 to state feed- 
back block 41 2 to update the detectiori state of the decoder. State feedback block 41 2 provides the new state as signal 
413 to decoder 409 for use in determining the detection thresholds for the next set of detection values 408. 

Rgure 8 illustrates a ML detector wherein the ML estimator comprises several FIR filters used to provide the linear 
. decision functions, Fx. Fy. etc. As shown, the equalized signal 405 is provided to the input of each FIR filter (800 and 

35 801 shown). The outputs of these filters, labeled Fx arid Fw, make up signal 408 which is provided to decoder 409. 
Within decoder 409, each value in signal 408 is compared with a state determined threshold. Further, a sign inversion 
is performed between positive phase and negative phase comparison operations. 

Each linear decision function in ML estimator 407 may be implemented as shown for FIR filter 800: A tapped delay 
line provides successive samples of the Input. Xo through y^: Each tap is multiplied by a coefficient value. Gq tiirough 

40 C4. respectively. Some of these coefficient values rriay be zero. After multiplication by the respective coefficients, the 
tapped signals are summed to provide the decision function estimate, is.g., Fx. An extra delay element and coefficient 
nrultipller C^i are shown to indicate the optibrial Inclusion of a look-back sample in the decision function of some 
embodiments. ' - ■ . - 

ML detector 406 consists of several function calculators embodied in ML estimator 407. and several comparators 

45 and a state machine embodied in decoder 409 and state feedback block 412. The transition rules between states of the 
state machine depend upon the logic results of cornparisons between the function values of the data samples and cer- 
tain state-dependent thresholds. Because there are fbuir non-zero samples in the equalized response of a single tran- 
sition, four bits are required to' represent the state of decoder 409. The state of the decoder is defined as the last four 
binary logic levels of the detected signal current, and is stored in state feedback block 41 2. Because of this format, state 

so feedback block 412 can be implemented with a four-bit shift register. The four bits can represent up to sixteen states. 

For a (1.7) RLL code EEPR4 channel, tiiere are ten legitimate NRZ (non-return-to-zero) states as shown in Figure 
3. wherein each transition indicates a transition from "6" to "1 ". or vice versa. However, due to the (1 ,7) RLL codingmin- 
imum for transition spacing, only six of the states require a decision to determine their next transition state and the 
decoded data. For example, state "0001" can only progress to state "0011." States "0010." "0100." "0101." "1010," 

55 "1011. "and "11 01 "are invalid NRZ states for (1.7) RLL coding. 

The task of ML decoder 409 is to determine whether a peak exists between clock cycle 1 and 2. which in turn deter- 
mines the next state. Figure 5 describes the decoding algorithm for tile detector. ML estimator 407 first computes the 
results of several linear functions (Fa. Fx, Fy and Fw) with the sample values yo, yi. ya. Ys and y^. and. in one embod- 
iment, look-back value yam. The results of these functions are compared with an equal number of state-dependent 
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thresholds in decoder :409^ The logic resuBs:of these comparisons are; assigned to the decision rvariables. A, X. Y jand 
W/Based on the values of ^:tbese:decision Yariables ard the cu/rent stateraoboCodo.^^ decoder determines. what. the 
next state and the: decoded, data are- according to the.rulesrin [Figure St^vT^ie. decisions for states 0001 ,1001 .1.1 10 and 
01 10 are not given in Figure 5 because their next state is^constrained by the (1 ,7)^coding scheme. 

5 Sample value y2 does not appear in equations [4] arid [Sl because it is cornmpri to both the signal and the minimum 

distance error waveform (as shown in. Figure 1). .TFhus. inclufiing the y2 sample would not iiicrease the detection thresh-, 
old, but would increase the tote^l addjtiye white Gaussian npi^ (AVVGNXby adding the ..noise of another sample. How- 
ever, the y2 sample do^ ^affect the noise, oprrelation and:.ttie in the ML 
estimator for the apparatus of the present invention improves the detector perfprrnance. , - ^ . ^ 

10 Further, a 'look-back" sample, (y^i) can be included as well to help. cancel noise This look-back sample 

can also, compensate for nonideal .equalizatlpri where, the actuab.sarnpie, value ris;^.n^^^ If y^i is 

included,. the four data state threshold feedback is optimally increased to f iye data stetes to set the proper comparator- 
thresholds. The trellis diagram of Figure 3 increases, from teri states ta sixteen states, which is .^ulyalent to increasing 
the partial response order to E^PR4. 

75 In Figure 5. function Fa is a baseline checking function. H determines whether the sarnpje values represent s signal 

that exceeds the minimum magnitude. threshold needed to indicate a peak caused by a transiti^ due 
to noise in the channel are thus ignored. Functiori Pa niust be true if a pe^ :r ^ . 

Fuctions Fx, Fy and Fw are responsible for determining whether the sample values fit a predeterrnined peak model. 
Fx performs a test for satisfaction of minimum conditions to distinguish a peak from the nearest valid sample model 

20 which does not form a peak. .Fy perfprms the same test as .Fx* but with a more stringent threshoid^^^ for Fy 

to be true, the sample model must be sufficiently close to an ideal peak model so as to differentiate it ifrpm any'other 
possible combination. Even a peak distorted by a subsequent peak will not pa^ 

Fw tests for sample combinations typical for a peak distorted by nearby peaks. Thus. for. an undistorted peak. Fx 
and Fy will be true, and for a peak distorted due to a nearby transition. Fx and Fw will be true Therefore, the resulting 

25 "peak detected" binary decision is A & X & (y + W).^ If this binary dedsipnjs.,true. then a trjansitiqn will be indicated at 
the output and in the state machine, unless such a transition is prohibited by the (1 ,7) coding. . 

The decoded data indicates whether there is a transition: at clock cycle zerOi that is. whether there is a transition 
from state bit c^ (note. Ci=do) to state bit dt. This decoded data is cprrpatible with most of the (1 .7) RLL decoders. The 
threshold values listed in the table of#igure 5 are considered nominal for an ideal Isolated EEPR4 pulse that is-"nor- 

30 malized" to have the sampled sequence value of 0, 1 .3.3. 1 ,0. Performance con be improved by fine-tuning these thresh- 
olds to give maximum margins in all condilipns.'^TTiis tuning, couldjeduce the . 

Vi denotes the sampled signal value at the ith dock cycle, and ajbjCjdj. represents the state of decoder 409, where 
aj. bj. Cj, dj are the last four binary logic levels of the data. In the following description, ye represents the current sample 
value that is being decoded, conresponding to .the clock cycle i=0. The symbols yi , y2. .Vs, and y4 are four look-ahead 

35 samples. The symbols aoboCodo represents the cun-enf state of the decoder and a^b^^ is the next state. From the- 
state value aobocodo. decoder 409 determines what the next state and, the decode<d dat are, according to the rules in 
Figure 6. A phase check function is not showyn;, !f a wrong polarity:Sig|nal: is detected the detector is reset to look for the 
proper phase. - ... >'r-%'...' •. ^ \ ■■■:y':rr^::.i .'. 

As previously described, a set of FIR filters maxirnizes the detection threshold between all possible signal wave- 

40 form patterns and their minimum distance en-or waveforms. The ykreighted sum outputs of the FIR f ilters are compared 
to a set of thresholds decided by the previously detected data tp^rnake the dedsipn fbr the cun-ent detection, the use 
of the previously detected data is analogous to the operation that pccurs in a decision feedback equalized detectoi- 
(DFE). In fact, the detector performance approaches that of fixed delay tree search (FDTS) with decision feedljack. The 
advantage of a PR channel over a: DFE is that it is easier and less, cqmplicated to^equalize to a PR target shaped pulse 

45 than to remove ISI. The disadvantage of PR is that its equalizatiprncprre^^^ - 

The FIR filter coefficient values effect the number of thresholds and hence the detector complexity. Coefficient val- 
ues also affect the complexity of a digital implementation. For example, rnultiplying by 0.5 is accomplished with.a digital , 
shift operation, and multiplying by .0:75 is aoconriplished with two shift operations and. an add operation. The number of 
equations may also affect the complexity. These constraints are. considered durii:»g the optimtzatbn of the coefficient.. 

50 values by. for instance, confining the alkjwable coefficient comt>inations under consideration to those meeting the con- 
straints. ■ - „.}"■■..■""- f:' ■ • • - \ ' [ ' . 

The optin^zation of the FIR filter coefficients of ML estimator 407 is illustrated in the flow diagram of Figure 9. The 
optimum coefficients are achieved by maximizing the, signal-to-noise ratio (SNR) at the output of the ML estimator FIR 
filters individually The signal value is, represented by the minimum distance error value for the respective FIR filter, 

55 whereas the noise value is represented by the square root of , the noise variance at the output of the FIR filter. Both of 
these elements, the signal value and the noise value, are functions of the FIR filter cojBfficients. Whereas the signal 
value may be calculated directly, determinatipn of. the noise wiance requir<^ euialysis of the PR chanriel frequency 
response, . -.■ 

In block 900 of Figure 9, the transfer furiction of; the PR. channel, including elements between the input off the 



7 



BNSDOCIO: <EP_0801484A2J_> 



EP0 801 484 A2 



10 



15 



Lorentzian pulse and the output of the target PR pulse, is determined. Once the transfer function of the PR channel is 
determined, in block 901. for the required user density, the noise autocorreilation coefficients are calculated from the 
equalized noise power density (i.e. the power spectrum of white noise filtered by the transfer function required to shape 
the Lorentzian pulse into a targeted PR pulse). In block 902. using either an exhaustive search or qualified selection 
process (e.g. a search limited to coefficient values requiring only shift and add functions etc.), a set of coefficients is 
selected for analysis, • ■ . .. . ^ 

In blocks 903 and 904, in any order, the noise variance at the output of the ML estimator and the minimum distance 
error value are calculated using the selected set of coeffidents. Next, an SNR gain value is determined from the mini- 
mum distance error value and the noise variance value in block 905, If, in subsequent block 906, sufficient coefficient 
sets have been analyzed, the set of coefficients providing the largest SNR gain value is selected in block 907. If further 
coefficient combinations remain to be analyzed, .the process, returns to block 902. 

The autocorrelation coefficients require that the variance of the noise power spectrum and the noise autocon^e- 
lation function a>(t) be calculated. The derivation of and <l>(t) follows. The Lorentzian pulse spectrum is given by the 
following equation: 



20 and the ideal equalized wave spectrum from the partial response equalizer Is given by: 

G(ij>)=TsO+D)" = Ts(1+-e^''^^)" 



[8] 



25 where <d is the frequency. PW50 is the pulse-width at fifty percent magnitude, and Ts is the sample period of the wave- 
form. . . 

Applying equations [7] and [8] to the additive white Gaussian noise at the input of the PR channel, the spectral den- 
sity of the colored noise at the output of the PR equalizer is: 



30 



\E{a>)\ 



2Ts 



^ 2 



50 



[9] 



35 The variance of the noise power density at the input of the ML estimator is derived from the above equation as fol- 
lows: - ■ >• -J ■•: -\ ^ \ ' • 



40 



'da> 



[10] 



where n/Ts is the bandwidth of G{a>)/E(a>). The variance integral can be reduced to the following summation: 



45 



2 4A/0 

KJ = — s 



[11] 



so where for n=0 (dicode), 6o=0.5; 
for n=1 (PR4). 6o=1 . 6i=1 : 

fbrn=2(EPR4), 61^3, 51=4,62=1 : and . . 

for n=53 (EEPR4). So=10. 5i=15. 82=6. 63=1, and 

55 where K = PWsq/Ts. For the (1 ,7) RLL code wherein the sample period or code period. Ts. is equal to two-thirds of the 
output NRZ data period (i.e. 2/3 bits per symbol). K is equal to 1 .5 times Du. the user density. 
The noise autocon-elatton function is given by the following inverse Fourier transform integral: 
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•2 75 1>^2^ 



50 



n?] 



which can be solved to obtain the following summation: 



10 



15 



2No 



[13] 



where T ^ =[(t/rs)+i] . and Tg =[(t/Ts)-i] . The autdcbrrelation coefficiente are obtained from equations [1 1] and [13] by 
dividing the autocorrelation function by the variance of the noise power density, p(t) = <]>(t)/cF2 . 

From the autocorrelation coefficients and the variance of the colored noise signal at the input of the ML detector, 
20 the noise variance at the output of each decision function Js calculated bajsed on the FIR coefficients of the respective 
decisions. The noise variance at the output of the FIR filter is given by: 



25 



30 



p n -i nn 

VAFA Y. o,nA = X Z Cmc,N,.CjNj] 

n n n n 



[14] 



n 



p n n /-1 * ■ - -j 

L/— 1 /-oy— 1 *; J 



where ^ is the noise variance prior to the FIR filter, p(t) is the autocorrelation coefficient at time "t", and Cj is the FIR 
filter coefficient for the i*^ sample. The FIR filter is of sample length "n". but the range of 1- begins at "-1 " to include the 
35 look-back sample (cM-tYM^) of one embodirnent. Wthout a look-back sample; the range lor "j" is between zero. arid "n". 
The signai-to-noise gain obtained by altering the FIR coefficient values can be measured by calculating the ratio of 
the minimum distance error, dMIN, dvided by the square root of the noise variance at the FIR filter output. 



40 



JmR 



[15] 



where 



45 



[16] 



so and Si and S2 are pulse signal waveforms for a minimum distance event. For relative comparison of SNR values, mul- 
tiplicative constants (e.g. cs^ in VAR) can be ignored. 

A larger value of dMIN signifies larger separation between minimum distance signal waveforms, providing more 
robust decisions. A smaller value of noise variance signifies a lower^probability that a false detection will occur due to 
noise for a given dMtN(FX). Therefore, by maximizing the SNR value above, performance of the ML detector is opti- 

55 mized. 

Because the number of coefficient. values is large for the determination of a maximum SNR value through deriva- 
tives, it is more straightforward to perform a cprnputer search for the optimal set of coefficients. A computer search pro- 
gram defining the range of coefficient values and the desired step sizes can evaluate the SNRs for a large number of 
coefficient combinations quickly and accurately. The best combination of coeffidents thus determined is then imple- 
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mented in the FIR filters of the ML estimator. Examples of programs for accomplishing this type of search are provided 
as Appendix 1 and Appendix 2. wherethe routine of Appendix 1 determines autocorrelation coefficients, and the routine 
of Appendix 2 determines the SNR values. The maixtrhufTl SNR value can be determined from the SNR values using a 
software sorting routine, or by reviewing the list of SNR values manually. 

5 For Du ranging from 1.5 to 3.0, a five-sample or six-sample jyiL look-ahead detector with .optimized coefficients 

approaches or outperforms Vrterbi detectors, especially in the higher usertJensity range. The detector performance can 
be optimized at a given Du. However, if a fixed set of ddefflciehfs is optimized for DU=2:0, the SNR performance is vir- 
tually the same as that for a Viterbi detector for Du ranging from 1 .6 to 3.0. If the look-ahead sanples are increased, 
detector performance will improve. Also, the detector linear equation coefficients can be optinriized to compensate for 

10 ndse coloring fronri sources other, than partisU resp^ 

Several variations off perfornnance versus ibbk-ahead ^rhple length and linear equation complexity are presented 
in Figure 7. Figure 7 shows detector signal-to-noise (S/N) gain versus Du when a Lorentzian pulse shape is ideally 
equalized to EEPR4. The gain of an ideal Viterbi Detector is shown ais a i-eference for the other detector an:angements. 
At Du=0. the noise is AWGN. 

15 Detector I is a five sample, three equation look-ahead detector optimizedlbr AWGN. Detector II adds the ya sample 
to compensate for correlated noise. y2 improves its performance by 0.56 dB at Du=2.0 and by 0.64 dB at Du-2.5. 
Detector ill is a more complex five sanrtple, three equation look-ahead detector whose performance is within 0.3 "dB of 
an ideal Viterbi detector for Du ranging from 1 .5 to 2.5. Detector IV is a five sarnple, three equation detector that utilizes 
the lookback sarrple and sixteen states. Detectbr V is a five sample; four equation look-ahead detector whose perform- 

20 ance exceeds the Viterbi detector over a wide range of Du. Detector Vl Is a . six sarhple, four equation detector with even 
greater potential performance improvement, - ^-i' ^ ; ^ 

A sample computer routine in BASIC for generating autocorrelation coefficients is provided as Appendix 1 . Further, 
a sample computer routine in BASIC lor generating a list of signal-tO:noise gain values is provided as Appendix 2. 
Thus, a method arid apparatus for reducing noise cbrrelation ih a partjat response channel has been described. 
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APPENDIX 1 



10 Tardal Response Noise Autocorrelation. Gpeffi 
20 PI=3,14159 
30 TR Selection 

40 -N^a :A(0)=.5 TDICODE: 
50 'N^l :A(0)=1 ^:A(1)=1 " TR4 

60 •N=2 :A(0)=3 :A(1)=4 :A(2)=1 'EPR* 
70 N=3 :A(OMO :A(1)=15 :A(2)=6 :A(3)=1 •EEPR4 
80 INPUT "K^";K :K2=K'^lC . 

90 B=EXP(K*PI) . ^ 

100 'Variance 
110 V=0 

120 FORI=0TON 

130 V=V+Aa)*2*(B*(K»c6sa*PI)+I^INa»PI))-K)/(I^ 
140 NEXT 
150 'Autocorrelation Coefficient 
160 INPUT Ts-jT 
170 P=0 

180FORI=0TON 
190 T1=T+I:T2=T-I 

200 P=P+Aa)*(B*K'^OS(Tl*PI)-K+Bni»SIN(Tl*PI))/(K2+Tl*Tl) 
210 P=P+Aa)*(B'^K*COSCr2'^PI)-K+B"^*SINCr2*PI))/(^ 
55 220 NEXT 

230 PRINT P/V :GOTO 160 . 
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■ 10 '1,7 £EPR4' ML s/N Gain'sear(ii;.;:: ; .f,: ' r'^- vVc!: ' . . ■ ' 

20 "Five Sample Look- Ahead Ghaninejl -With Golpre . ^ - . - 

30 Ideal EQ Filter- - - ^ ' . r . n .; . : ,. r ".. . - . 
40 •(Y0,Yl,Y2,V3;Y4)=(A,B;CaE) - • 

50 *SignaI=(U^4,-l)-(04^,2,-2)=a^AriA)jS^^*A+4 . 

60 FORA=.dT6i ' -■ '-v 

70 FOR B=.5 TO 1.5 STEP .1 ci-.-: :-:^-^r: ' \ 

80 FOR C=-.5 TO 1.5 STEP.l 

90 FOR D=-l TO .0 STEP .1 = - - ^ ^ : : • - ^ 

100 FORE=.5TOi.5 STEP.l .. .<:■■■■■- r- -^^ -yy- 

110 s=A»A+4»B*B+b**p+E»E rJ /■-:^V:V^\.J i ■;•■■: ...■::„.':;.'". 

120 Rl=.4826 :R2=-.1722 :R3=>ia . . ' : 

130 V=A»A+B*B+C*G+D»D+E'^E+2»((A»B+B»G+,G'tD+Dt"*E)*Ri .^ .s 

+(A*G+B»D+C»E)*R2+(A»D+B*E)»R3+A»E'^R4) 
140 SN1=SQR(S)/SQR(V) : . v r ? ;.- 

150 Rl=.2999 :R2=-.3907 :R3=-.i524 :k4=.iD324 •I>u=1.5 

160 V=A»A+B»B-»<:»C+D»D+E»E+2»((A»B+B»C+C*b+t)*E)*Rl 

+(A»C+B»D+C*E)»R2+(A»D+B»E)»R3+A»E»R4) 
170 SN2=SQR(S)/SQR(V) ; v 

180 Rl=.1077 Jl2=-.5348 :R3=-.021 :R4=.086 •Du=2 
190 V=A»A+B»B+C»C+D*D+E»E+2*((A»B+B»C+C*D+D»E)*R1 
+(A»C+B»D+C»E)»R2+(A»D+B»E)»R3+A»E*k*) 

200 SN3=SQR(S)/SQR(V) 

210 Rl=-.0726 :R2=-5852 :R3=.1631 :R4=.0787 '0^=2.5 
220 V=A»A+B*B+C'^+D»D+E*E+2*((A»B+B*C+C*D+D*E)»R1 

+(A*C+B*b+C*E)*R2+(A»D+B*E)*R3+A*E*R4) : 
230 SN4=SQR(S)/SQR(V) T ' 

240 Rl=-.2286 :R2=-.5547 :R3=.3388 d^W.Q003 'Du^S. . 
250 V=A*A+B»B+C»C+D»D+E»E+2^(A»B+B»c4c*r>+f)»E)»Rl . 

+(A»G+B*D+C»E)»R2+(A»D+B»E)*R3+A*E»R4) 
260 SN5=SQR(S)/SQR(V) 

270 PRINT USING "####.##"^A;B;G;D;E;20»LOG(SN1)/LOG(10); . 

20»LOG(SN2)/LCK3(10>;20*LOG(SN3)/LOG(10);20*LOG(SN4)7lX)G(lp); 

20»LOG(SN5)/LOG(10) ~ ~ 

280 NEXT:NEXT:NEXTa^^EXTJ^^EX^ . •. . . - 



Claims 

1. A partial response channel comprising: 

an equalizer having a first frequency response for shaping an input pulse signal into a target pulse signal, said 
first frequency response coloring an adcfitive noise component of said infiut pulse' signal and generating a cor- 
related noise component of said target pulse signal; 

a detector receiving said target pulse signal and providing a binary data signal, said d^ector comprising a plu- 
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rality of detection decision circuits for receiving said target pulse signal and providing a plurality of decision val- 
ues, wherein at least one of said plurality of detection decision circuits comprises a finite impulse response 
(FIR) filter having a set of coefficients maximizing a slgnal-to-noise ratio; 

a decoder having a plurality of state dependent threshold values, said decoder providing said binary data sig- 
5 nal based on comparisons between said plurality of decision values and said state dependent threshold val- 

ues. 

2. The apparatus of claim 1 wherein said signal to noise ratio comprises a minimum distance error value of said target 
pulse signal divided by a square root of a noise variance at an; output of said FIR filter. 

10 

3. The apparatus of Claim 1 or Claim 2, wherein said FIR f ilter comprises a tapped delay line providing a plurality of 
tapped signals, each of said tapped signals being multiplied by a coefficient value from said set of coefficients to 
provide weighted signals, said weighted signals being summed to provide one of said decision values. 

15 4. The apparatus of any preceding claim,. wherein said equaliser has a first frequency response that is a lowpass func- 
tion with high frequency boost. * * 

5. The apparatus of any pi-ecedihg claim, wherein said diecoder further comprises a state machine wherein a transi- 
tion between a first state and a secorid state is dependent upon said cornparisons. 

20 

6. The apparatus of any preceding claim, wherein the decoder provides a binary data signal that is run-length-limited 
(RLL) encoded. 

7. The apparatus of any preceding claim, wherein said FIR filter has a second frequency response arranged for whit- 
es ening said correlated noise component of said target pulse signal. 

8. The apparatus of any preceding claim, wherein said detector is a maximum-likelihood (ML) detector. 

9. A method Ibr reducing noise correlation in a partial response (PR) channel comprising a noise correlating element 
30 coupled to a detector, said detector further comprising a finite impulse response (FIR) filter coupled to a decoder, 

said method comprising the steps of: . f " 

determining noise autocorrelation Coefficients for a correlated noise signal at an input of said FIR filter; 
determining a plurality of noise variance values at an output of said FIR filter for a plurality of FIR coefficient 
35 combinations; 

determining a plurality of minimum distance error values of said FIR filter for said plurality of FIR coefficient 
combinations; ... 

utilizing an optimal coefficient combination from said plurality of FIR coefficient combinations which maximizes 
a signal-to-noise ratio (SNR), said SNR comprising a function of said minimum distance error value and said 
40 noise variance value for each of said plurality of FIR coefficient combinations. 

10. The method of Claim 9 further conrprising; determining said SNR for each FIR coefficient combination by dividing 
an associated minimum distance error value by a square root of said noise variance value. 

45 11. The method of Claim 9. whereiasaid step of determining said noise autocorrelation coefficients further comprises: 

determining a variance of a noise power density of said correlated noise signal based on a transfer function of 
said noise correlating element; 

dividing an autocorrelation function by said variance of said noise power density. 

so 

1 2. The method of any of Claims 9 to 11 , wherein said step of determing said noise variance values comprises deter- 
mining said noise variance values for a selected user density (Du). 

13. The method of any of Claims 9 to 12. further comprising limiting said plurality of FIR coefficient combinations to 
55 coefficients requiring only shifts and adding operations. 
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Positive Phase 


Funcdori of Look-Ahead 
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Nonunal Thresholds at State 
0000 _ 1000 ^ 1100 


Decision 
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Baseline Check 

Fa = yi+y2 


Fa>23 


Fa>2^ 


Fa > 1.5 


A 


Peak Check 

Fx = o^yo+yi-ya 

Fy = 03yo+yi-y3 
Fw = 03yo+yi-H)5y2-0^y3+y4 
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. . . Colorfed Noisfe S/N Cain" " 


'User Density,-Du 


Detector Tvoi» 


p 
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1.5 


2J0 


2.5 


3.0 
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10 


7.76 
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8.20 
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1. 5 Sample, 3 EquaLtiqn 
Fx«Fy«05yo*rt-y3' , 




7,75 i 


7*40 


7.13 


7.04 


7.43 


5 Sample, 3 Equation 
Fx « Fy a 05yo+yi-y3 \ 
Fw = OJ5yo+yi +0 Jy2r03y3+y4 




7J9 . 


7.64 


7j69 
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Fx = Fy = 0.8yo40.6yi+03yary5 
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^ IV. 5 Sample, 3 Equation 
Fx = Fy = 0.4yMl+0^4yo40-7yi+0.2yH}^y3 
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V. 5 Samplj^ 4 Equation 
Fx - Fy » 0.75yo+yi-y3-075y4 
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Fv « yD-fr03yi+0.75y2-y3+0-5y4 

i I 
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VL 6 Sample, 4 Equation 
Fx - Fy « Fz « 0-75yo+yi-y3^^75y4 

: Fw«yo+0.75yi+y2-y3+y4 
Fv = 0-75yo+yi+0 Jyi^^y340^y4+y5 
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9^ 


10.90* 



1. Optimized for Du»2.0 

2. Optimized at each Du 
- (Fa « yl+y2, for all above cases) 
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(54) Improvements In or relating to signal processing 



(57) A method and apparatus for reducing noise 
correlation in a partial response channel through optimi- 
zation of a look-ahead maximum likelihood (ML) detec- • 
tor. In the method of the present invention, the ML 
detector (406) is optimized in light of the noise correla- 
tion generated by the partial response channel. The 
improved ML detector provides comparable perform- 
ance to. or better performance than, a Viterbi detector in 
the presence of colored noise. In the present invention, 
a set of finite impulse response (FIR) transversal filters 
(800.801) are used as the ML est'mator(407) for the 
look-ahead detector. The weighted sum outputs of the 
FIR filters are compared to a set of thresholds based on 
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previously detected data to make the decision for cur- 
rent detection. The present invention improves the ML 
detector's performance and reduces its complexity by 
optimizing the coefficients of the FIR filters in the pres- 
ence of the correlated or colored noise. The SNR (sig- 
naf-to-noise ratio) of each FIR filter is determined for a 
range of coefficients based on the noise autocorrelation 
of the channel for a given user density, and the coeffi- 
cients providing the highest SNR are selected for each 
decision function. The result is a noise whitening ML 
detector providing innproved performance and lower 
complexity than prior art ML detectors. 
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